An unusual stellar death on Christmas Day 
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Long Gamma-Ray Bursts (GRBs) are the most dramatic examples of massive stellar deaths, usually as- 
sociated with supernovae LU- They release ultra-relativistic jets producing non-thermal emission through 
synchrotron radiation as they interact with the surrounding medium f2l. Here we report observations of 
the peculiar GRB 101225A (the "Christmas burst"). Its 7-ray emission was exceptionally long and followed 
by a bright X-ray transient with a hot thermal component and an unusual optical counterpart. During the 
first 10 days, the optical emission evolved as an expanding, cooling blackbody after which an additional 
component, consistent with a faint supernova, emerged. We determine its distance to 1.6 Gpc by fitting the 
spectral-energy distribution and light curve of the optical emission with a GRB-supernova template. Deep 
optical observations may have revealed a faint, unresolved host galaxy. Our proposed progenitor is a helium 
star-neutron star merger that underwent a common envelope phase expelling its hydrogen envelope. The 
resulting explosion created a GRB-like jet which gets thermalized by interacting with the dense, previously 
ejected material and thus creating the observed black-body, until finally the emission from the supernova 
dominated. An alternative explanation is a minor body falling onto a neutron star in the Galaxy Q. 

On Dec. 25, 2010, 18:37:45 UT (Tq), the Swift Burst Alert Telescope (BAT, 15-350 keV) detected GRB 101225 A, 
one of the longest GRBs ever observed by Swift IH with T90 > 2000 s (the time in which 90% of the 7-ray energy 
is released, O). A bright X-ray afterglow was detected for two days and a counterpart in ultraviolet, optical and 
infrared bands could be observed from 0.38 h to two months after the event (see the Supplementary Information, 
SI). No counterpart was detected at radio frequencies ||6l|3. 

The most surprising feature of GRB 101225 A is the spectral energy distribution (SED) of its afterglow. 
The X-ray SED is best modeled with a combination of an absorbed power-law and a black-body (BB). The 
UV/optical/NIR (UVOIR) SED (see Fig. 1) can be modeled with a cooling and expanding BB model until 10 
days post burst, after which we observe an additional spectral component accompanied by a flattening of the light 
curve (Fig. 2). This behaviour differs from a normal GRB where the SED follows a power-law due to synchrotron 
emission created in shocks when the jet hits the interstellar medium (e.g., II2II). 

An optical spectrum taken two nights after the burst does not show any spectral lines (see SI). Were fit the SED 
and fight curve with the template of SN 1998bw, a SN Type Ic associated with GRB 980425 |44|, and obtain a 
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redshift ofz = 0.33. At this distance, the SN has an absolute peak magnitude of only My = —16.7 mag, which 
makes it the faintest SN associated with a long GRB ||9l[T0]|. In contrast, the 7-ray isotropic-equivalent energy 
release at z = 0.33 is > 1.4 x 10^^ erg, typical for other long GRBs but more luminous than most other low- 
redshift GRBs associated with SNe 1 1 1 1. We detect a possible host galaxy in g' and r' bands with OSIRIS/GTC at 
6 months after the burst with an absolute magnitude of only M^^abs = —13.75 mag, ~ 2 mag fainter than any other 
GRB host |[T2ll . Although its blue colour matches that of a star-forming galaxy, our observations do not allow us 
to resolve it as an extended source. 

Atz = 0.33, the X-ray BB has a radius of - 2 x 10^^ cm (-1 R©) and a temperature of - 1 keV (10^ K) 
at 0.07 d with little temporal evolution. Such a thermal component, attributed to the shock breakout from the star, 
had also been observed for XRF 060218 d, XRF 100316D/SN 2010dh (H and GRB 090618 (HI, all nearby 
GRBs associated to Type Ic SNe (lUITTlllSinil, with similar temperatures but larger radii. The UVOIR BB starts 
with a radius of 2 x 10^^ cm (~13 AU) and a temperature of 8.5 x 10^ K at similar times and evolves considerably 
over the next 10 days reaching a radius of 7x lO^^cm and temperature of 5,000 K. The evolution of the two BB 
components suggests that they must stem from different processes and regions. 

An appealing model is a Helium - neutron star (NS) merger with a common envelope (CE) phase, a model 
that had been proposed earlier as a possible progenitor for GRBs |[l9l |20| [211 . In this scenario, a binary system 
consisting of two massive stars survives the collapse of the more massive component to a neutron star (NS). When 
the second star leaves the main sequence and expands, it engulfs the NS, leading to a CE phase and the ejection of 
the hydrogen envelope and part of the helium core as the remnant spirals into the center of the second star. When 
the NS reaches the center, angular momentum forms a disk around the remnant of the merger, allowing for the 
formation of a GRB -like jet. The remnant might be a magnetar whose prolonged activity can explain the very long 
duration of the GRB. 

The interaction of this ultra-relativistic, well-collimated jet with the previously ejected CE material can explain 
both the X-ray and UVOIR emission components. Estimating that the inspiral takes 5 orbits or 1.5 yr and material 
is ejected at escape velocity, the outer ejecta are at a few 10^^ cm at the time of the merger, consistent with the 
radius of the UVOIR BB. We assume that the ejecta form a broad torus with a narrow, low-density funnel along the 
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rotation axis of the system that permits the passage of the 7-radiation generated in the jet. Most of the jet hits the 
inner boundary of the CE-ejecta and only a small fraction of it propagates through the funnel. The X-ray emission 
is produced by shocks created by the interaction of the jet with the inner boundary of the CE shell. As the jet 
passes through the funnel, it decelerates due to the increased baryon load and shear with the funnel walls so that 
no regular afterglow signature is produced. When the now mildly relativistic, mass-loaded jet breaks out of the 
CE-ejecta, it produces the UVOIR emission in the first 10 d. As the supernova shock expands beyond the CE shell 
we observe a small bump in the lightcurve at ~30d. The He-NS merger scenario naturally assumes a relatively 
small Ni-production, leading to a weak SN. 

A somewhat similar scenario might also explain XRF 060218 albeit with a different progenitor system pro- 
ducing a brighter SN and a fainter GRB. The remnant might be a magnetar whose prolonged activity can explain 
the very long duration of the GRB. A thermal component in the optical can also be fitted for the X-ray out- 
burst associated with SN 2008D, a Type lb SN in NGC 2770 |22|, which, however, had no thermal component 
in X-rays. And last, it is possible that a GRBs shows a thermal component in X-rays (e.g., GRB 090618) but 
shows a classical, bright power-law shaped afterglow 1 18|. GRB 101225 A might be a member of a newly defined 
class of "blackbody-dominated", SN-associated long-duration GRBs arising in very dense environments created 
by the progenitor systems themselves, which thermalizes the high-energy output from the collapsing star. The 
non-relativistic, uncollimated emission from this scenario makes it difficult to detect them at higher redshifts. This 
makes GRB 101 225 A a fortunate coincidence that allows us to derive conclusions about the progenitor system and 
its environment by a new variety of massive stellar death, which had so far been only proposed to exist theoretically. 
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Figure 1: Temporal evolution of the UVOIR spectral energy distribution. The UV, optical and IR counterpart were detected 
by UVOT, the UV telescope onboard Swift and several ground based facilities from 0.38h to nearly 2 month after the GRB. 
This plot shows the evolution of the SED from the onset of the optical observations at 0.07 days to 40 days for all epochs with 
sufficient data to model the SED shape. Filled circles are detections, triangles mark upper limits. The orange line on top of 
the BB model at 2.0 days shows our flux-calibrated spectrum taken with the OSIRIS/GTC. The SED evolution requires two 
different components, a simple expanding and cooling BB up to 10 days and an additional SN component for the last 3 epochs. 
The solid line shows the combined evolution taking both the BB and the SN beyond 5 days into account, whereas the dashed 
line shows the evolution of the BB component alone. The UVOIR BB evolves from an initial temperature of 43,000 K (0.07 d) 
to 5,000 K (18 d) and increases in radius from 2x 10^^ cm to 7x 10^^ cm at the same timescale. The SN component is best fit 
with a template of the broad-line Type Ic SN 1998bw which was associated to GRB 980425. For the fit of the SN component 
we used the SED at 40 days. Reanalyzing UVOIR data of XRF 060218 [13] and SN 2008D |22 | we find a similar thermal 
component over the first 3-4 days after which the SN starts to dominate. The UVOIR BB of XRT 060218 shows a similar 
radius, temperature and evolution as GRB 101 225 A and is equally inconsistent with radius and temperature from the X-ray BB. 
SN 2008D shows a steeper evolution in the radius, consistent with the continuation from the shock breakout, a thermal X-ray 
component was not observed (see SI). 
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Figure 2: Light curves of GRB 101225A in X-ray and UV/optical/IR bands, (a) X-ray lightcurve, (b) light curves in UV 
to IR bands. Filled points are detections, triangles upper limits. The solid lines denote the combined light curve from the BB 
and the SN component, excluding the contribution from the host galaxy. The evolution of the BB component is shown as a 
dashed line, the SN as a dotted line. Observations started almost simultaneously in X-rays and optical/UV wavelengths. The 
X-rays reached a peak flux of 4.34 x 10~^ erg cm~^ s~^. After an initial shallow decay of slope ^-i iosio.oii 21 ks, 
the X-rays show a strong decay with a slope of ^-5-95±o.2o^ inconsistent with synchrotron emission. The UVOIR light curves 



have a shallow maximum at the beginning, with different peak times for the individual bands due to the maximum of the BB 
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emission passing through the spectrum. The second component emerging at around 10 days post-burst is the contribution of 
an underlying SN, modeled with the GRB-SN 1998bw as a template, stretched in time by a factor of 1.25 and decreased in 
luminosity by a factor of 12 (in restframe). The absolute luminosity of the SN is My = —16.7 mag, the faintest SN associated 
with a GRB. At ~ 180 days we detect the very faint host at magnitudes of g' = 27.36 =b 0.27 and r' = 26.90 =b 0.14 or 
Mabs = - 13.7 (0.001 L*). (c) Color image of the field of GRB 101225A observed at 40 days (indicated by a grey bar in panel 
b) with the afterglow marked by a box. 
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Supplementary information 



1 BAT data analysis and fitting 

GRB 101225A was detected by BAT on-board the Swift satellite 123 on Dec. 25, 2010 at Tq = 18:37:45 UT as 
an image trigger |24|. It was already in progress both when the source entered the BAT field of view and when it 
left the field of view due to Swift orbit-constrained slews |[25]|26l. Therefore, we can only give lower limits on the 
total burst fluence and the Tgo duration. The total fluence of the intervals covered in the observations adds up to 
(5.6 ± 0.7) X 10~^ ergcm"^ (implying a total energy release in 7-rays of ^j^o > 1.4 x 10^^ erg at z = 0.33.), 
which is a lower limit to the total gamma-ray emission. No emission was detected in a previous observation of the 
field at T — To = —4950 s and we therefore put a lower limit on the duration of T90 > 2000 s. This is one of 
the highest durations ever observed for a GRB, comparable to the longest burst observed by Swift, GRB 0904 17B 
(271. The BAT light curve is shown in Fig. [l] The BAT-observed peak flux of (3.25±0.47) x 10"^ erg cm-^ 
s"^ in the 15 - 150 keV range occurred in the interval T - Tq = +1372 to T - Tq = +1672 s. No other 7- 
ray instrument detected GRB 101 225 A, although the MAXI instrument on board the ISS (2 - 10 keV) reported a 
marginal detection at T — Tq = +1002 s coincident with the BAT position 1281 . 

The time-averaged spectra from Tq to T - Tq = +963 s and from T - Tq = +1372 to T - Tq = +1672 s 
are best fit by simple power-law models with photon indices ofF = 1.9 + 0.4 and 1.9 + 0.2, respectively. The 
total fluences in each time period in the 15 — 150 keV band taking these fits are (1.7 + 0.4)x 10~^ and (9.0 + 0.2) 
xlO~^ erg cm~^. All quoted errors are at the 90% confidence level. The BAT spectra are almost equally well 
parametrised by models using a cutoff power-law or a blackbody fit due to the low signal-to-noise ratio of the 
event. Epeak using a cutoff model is poorly constrained to 38 + 20 keV. The blackbody fit gives a temperature of 
kT = 10. 1 + 1.1 keV. Errors are at the 68% confidence level. 

We also examined the BAT data to search for persistent emission after the trigger. For this we used the daily 
sky image mosaics produced as part of the BAT hard X-ray transient monitor which cover a single energy band 
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Figure 1: BAT 7-ray lightcurve in the 15-150 keV band. The sum of the flux was obtained using a fixed power- 
law index of F = 1.87 from a fit to the most intense part of the burst. The burst started before the beginning of the 
BAT data at ~ T — Tq = — 100 s and probably continued while the source was not in the BAT field of view from 
T - To = +1091 to T - To = +1372 s. The latest upper limit before the burst was 2.65 x 10"^ erg cm-^ s"^ at 
T - To = -4950 s. Error bars are at 90% confidence. 

of 15 - 50 keV. We found a 5.3cr excess (0.0048 ± 0.0009 count cm-^ s-^) on Dec. 25, 2010 (MJD 55555), the 
day of the trigger, and a positive excess in the count rate (> la or 0.001 1 count cm~^ s~^) over the next ten days 
(until MJD 55565). We determine the probability that such a sequence of excess rates would occur by chance. To 
do this, we examine the light curves of 106 "blank sky" points tracked in the BAT transient monitor. These are 
points chosen randomly across the sky at least 10 arcmin from any known X-ray source. Any positive flux from 
these points is expected to be due to chance fluctuations. In these 106 light curves (> 200, 000 data points), we 
find only one sequence of six consecutive days showing a positive excess and none with more than six days. The 
chance probability of ten days of excess flux is less than 1/200,000, so the observed prolonged emission is likely 
real. 

2 XRT data analysis and fitting 

The Swift XRT data were processed with version 3.7 of the XRT data reduction software (released as part of 
HEASoft 6.10 on 2010-09-28) and the corresponding calibration files used for subsequent spectral analysis. 
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The object was detected by XRT from 1.4 ks to 10^ s after the trigger. Data were collected in Windowed Timing 
(WT) mode for the first 7.3 ks after the trigger followed by Photon Counting (PC) mode for the rest of the obser- 
vations. The peak flux in X-rays is 4.3 x 10~^ erg cm~^ s~^, the total observed fluence 8.2 x 10~^ erg cm~^ and 
the unabsorbed fluence 1.1 x 10~^ erg cm~^. At z = 0.33 this corresponds to a total energy release in X-rays of 
3.6 X 10^^ erg. Spectra were extracted for individual snapshots of data (one snapshot corresponds roughly to one 
orbit constrained by the observability of the object during the orbit) and were further timesliced into 100 s bins for 
the initial snapshot (1.4 — 1.8 ks after the trigger). 




Figure 2: Fit to the X-ray spectrum. Panel a): X-ray spectrum from XRT during the first snapshot, the dashed 
line indicates the contribution of the power-law component, the dotted line shows the BB component. Panel b): 
Ratio between the observed data and the fitted model. 

We tried a variety of fits to the X-ray data, using XSPEC version 12.6.0, with the result that an absorbed 
power-law plus BB component provided a good fit to the data (see Fig. [2]). For the fit we used the Tubingen-Boulder 
absorption model with the Wilms abundances 1291 and Verner absorption cross-section ISOl . As shown in Fig. [3j 
there is little spectral evolution within the first snapshot of data. The best fit taking the entire first 367 s of data is 
a power-law with a photon index F = 1.83to;lo' ^ blackbody of temperature 0.96 ± 0.13 keV (1.11 x 10^ K) and 
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Figure 3: Results from the fits to the first snapshot of the X-ray data. The panels show from top to bottom: a) 
count rate during the first snapshot, b) photon index F, c) BB Temperature in keV, d) the contribution of the BB 
to the total emission in percent, e) radius of the emitting BB at 2: = 0.33 and f) total absorbing column density in 
X-rays (the Galactic column density in the line-of-sight is 7.9 x 10^^ cm~^). 

a total absorbing column of (2.2 ± 0.3) x 10^^ cm~^, for a of 420.7 for 379 degrees of freedom. The Galactic 
column density in this direction is 7.9 x 10^^ cm~^. 

The inclusion of the BB is significant at the > 99.9999% level, the contribution of the BB to the total emission 
is around 20% (see Fig. [3]). The second snapshot of data (also in WT mode) is again better fit with a BB in addition 
to the power-law, with F = 2.18t^;J^, black-body kT = OM^olr keV and NH = (2.7 ± 0.2) x 10^^ cm-^, 
with x^/d.o.f = 378/421. This BB is significant at 99.987%. For the X-ray data after the second snapshot, no 
BB component is required and a simple absorbed power-law provides an acceptable fit, likely due to the lower 
signal-to-noise ratio at later times. 
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Table 1: Result of different model fits to the X-ray data of the first snapshot. PL is a pure power-law model, 
PL+BB a combination of a simple PL and a BB with one temperature while a PL+diskBB includes emission from 
BBs with different temperatures, PL+compt includes a comptonized component in addition to the power-law. F is 
the photon index of the power-law, kT the BB temperature, the column density is the total density, including the 
Galactic absorption. The last column shows the of the different fits and the F-test value compared to the simple 
absorbed PL model. 



Model 


F 


kT (keV) 


Nh (10^^ cm-^) 


X^/d.o.f. (F-test) 




PL 

PL+BB 

PL+diskBB 

PL+compt 


1.72±0.03 

1 OQ+0.13 

^•^■^-0.10 
1 79+0.36 

1.79±0.04 


0.96±0.13 
1.64±0.35 

<22 


0.24±0.02 
0.22±0.03 
0.18±0.04 

rvo.+0.03 

^•■^^-0.02 


468/381 

421/379(1.95x10" 
417/379 (3.19x10- 
424/377 (1.55x10" 


-9) 
-10) 



The detection of N// in excess of the Galactic column density can be used to constrain the redshift since 
column density and redshift are coupled. We fitted the spectrum of the first orbit fixing the Galactic absorption 
to 7.93 X 10^^ cm~^. In addition to the absorbed powerlaw plus BB component as described above we allowed 
for a third component leaving both the column density and the redshift free to vary (in contrast to the fit described 
above where the redshift was fixed to 0.3). The resulting contours are plotted in Fig. |4] The 99% upper limit for 
the redshift is 0.5, the 90% upper limit on the redshift is 0.35, consistent with our findings from the UVOIR SED 
fitting. 

We also checked for a possible periodicity in the X-ray data. To this end, light curves were extracted with 18 ms 
bins which is the best time resolution available for WT mode. Using the Kronos powspec tool, no significant 
periodic signal was identified with a frequency between 0.005 and 28 Hz (0.04 and 200 s) in either the first or 
second snapshot of data. The 3 a limit on the variation in flux of any periodic signal during the first WT snapshot 
is 6.3%. 

3 Uy, optical and IR data analysis 

3.1 Swiftl\]\OT 

Swift/UYOT began observing GRB 101225A 1373 s after the BAT trigger, simultaneous with the XRT observa- 
tions. The automatic target sequence did not start until the end of the BAT image trigger at 23 minutes. The 
optical counterpart was found to be blue, with strong detections in the UV filters (uvwl, uvm2, uvw2), weak 
detections in the b and u filters, and no detection in the v filter. The data were processed using the standard Swift 
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CO 




Nh (1022 cm-2) 

Figure 4: Redshift vs. absorbing column density confidence contours. The fit is based on the data from the first 
snapshot. Plotted are the 1, 2 and 3 cr contours together with the 90% and 99% confidence interval. 



software tool uvotmaghist within HEAsof t 6 . 9 and the latest calibration files (20101231). 

We extracted counts using a circular aperture with a radius of 5 arcsec for count rates above 0.5 counts s~^, 
and a 3 arcsec aperture where the count rate drops below 0.5 counts s~^, as well as a source-free background 
region. The tool uvotmaghist applies coincidence-loss corrections and aperture corrections. The count rates 
were converted to flux density using standard photometric calibration 131113211 . 

3.2 McD 2.1m 

The CQUEAN instrument (Camera for QUasars in the EArly uNiverse; Park et al. 2011, in preparation) on the 
2.1m Otto-Struve telescope at McDonald Observatory, Texas, USA, observed the optical counterpart starting at 
01:16:23 UT, on Dec. 26, 2010 (6.64 h after the burst). Three exposures of 300 s each were taken in r^ ^^ z' , 
and Y bands under photometric conditions. The data were reduced with standard procedures of dark and flat-field 
corrections. The afterglow is detected in the r^ ^^ and z'-band images, the Y band gives an upper limit only. 
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Table 2: Log of observations in the UVOIR range. Values are not corrected for Galactic extinction. 



Mid t-tQ 


Exposure 


Filter 


Telescope 






(days) 


(s) 










Premaging 


3 X500 


„/ 

8 


3.5mCFHT 


>26.9 (27.2ih0.5) 


<0.06 (0.048 ±0.22) 


Premaging 


3 X500 




3.5mCFHT 


>25.5 


<0.22 


0.01848 


168 


w2 


UVOT 


>21.36 


< 10.38 


0.07041 


1431 


w2 


UVOT 


21.56± 0.20 


8.65± 1.80 


0.17373 


6719 


w2 


UVOT 


21.63± 0.11 


8.06± 0.88 


0.30739 


6679 


w2 


UVOT 


21.76d= 0.12 


7.19± 0.85 


0.45280 


5805 


w2 


UVOT 


21.96± 0.15 


5.96± 0.91 


0.81302 


12039 


w2 


UVOT 


22.57± 0.17 


3.42± 0.58 


1.00869 


11753 


w2 


UVOT 


22.37± 0.16 


4.08 ± 0.67 


1.41736 


23440 


w2 


UVOT 


22.6 1± 0.20 


3.27 ± 0.66 


1.75211 


23368 


w2 


UVOT 


23.45± 0.30 


1.51 ± 0.49 


2.44862 


74516 


w2 


UVOT 


>23.73 


< 1.17 


4.07964 


138747 


w2 


UVOT 


> 24.20 


< 0.76 


7.52954 


377712 


w2 


UVOT 


>25.39 


< 0.25 


0.01818 


319 


ml 


UVOT 


> 20.81 


< 17.14 


0.07515 


1431 


m2 


UVOT 


21.97± 0.31 


5.89 ± 1.94 


0.61452 


899 


m2 


UVOT 


21.90± 0.21 


6.34± 1.33 


0.95369 


12104 


m2 


UVOT 


22.00± 0.15 


5.74± 0.83 


1.18487 


18396 


m2 


UVOT 


22.47± 0.23 


3.74± 0.87 


1.48860 


23468 


m2 


UVOT 


22.46± 0.19 


3.76± 0.73 


1.85629 


29528 


m2 


UVOT 


22.97 ± 0.22 


2.35 ± 0.53 


2.51507 


40973 


m2 


UVOT 


23.34± 0.30 


1.68± 0.53 


4.08285 


138701 


m2 


UVOT 


> 24.25 


< 0.73 


0.01846 


318 


wl 


UVOT 


>21.15 


< 12.64 


0.07752 


1431 


wl 


UVOT 


>22.10 


< 5.26 


0.65205 


5571 


wl 


UVOT 


21.81± 0.17 


6.88± 1.15 


0.96984 


16520 


wl 


UVOT 


21.72± 0.16 


7.46± 1.19 


1.37145 


18904 


wl 


UVOT 


22.23 ± 0.26 


4.65 ± 1.25 


1.71425 


29052 


wl 


UVOT 


22.46± 0.25 


3.76± 0.97 


2.44380 


74329 


wl 


UVOT 


> 22.97 


< 2.36 


4.07673 


138760 


wl 


UVOT 


>23.22 


< 1.86 


0.01789 


169 


u 


UVOT 


> 20.46 


< 23.68 


07228 




u 




zi.jyit u.zo 


8 42 ± 2 47 




103571 




UVOT 


22.33 lb 0.26 


4.24 ± 1.14 


2.44533 


74215 


u 


UVOT 


>21.82 


< 6.80 


4.07781 


138647 


u 


UVOT 


> 22.06 


< 5.46 


0.01817 


169 


b 


UVOT 


> 19.94 


< 38.37 


0.06802 


1430 


b 


UVOT 


>20.86 


< 16.51 


0.16319 


6726 


b 


UVOT 


21.53± 0.30 


8.83 ± 2.86 


0.30647 


8349 


b 


UVOT 


>21.83 


< 6.76 


1.14624 


127458 


h 


UVOT 


>22.19 


< 4.85 


2.44615 


74224 


h 


UVOT 


> 2 1.00 


< 14.47 


4.07834 


138632 


h 


UVOT 


>21.33 


< 10.64 


0.01789 


318 




UVOT 


> 19.35 


< 66.30 


0.07278 


1431 




UVOT 


>20.36 


< 26.12 














0.31290 


5819 




UVOT 


>20 69 


< 19 20 


1.11739 


21x180 


V 


1.23mCAHA 


22.47 ±0.19 


3.73 ±0.65 


1.18728 


121098 




UVOT 


>21.08 


< 13.39 


2.45087 


74279 


V 


UVOT 


>20.50 


< 22.93 


4.08128 


138537 


V 


UVOT 


> 20.92 


< 15.56 


39.11207 


6X180 


g' 


OSIRIS/10.4mGTC 


> 26.3 


< 0.11 


39.49403 


5X180 


g' 


GMOS/8mGemini 


26.80±0.35 


0.07±0.03 


--180 


42 X 200 


g' 


OSlRlS/10.4mGTC 


27.21 ±0.27 


0.047±0.010 


1.04545 


19X 180 




1.23mCAHA 


22.61 ±0.16 


3. 28 ±0.48 


0.29887 


3X300 




CQUEAN/2.1mMcD 


22.43 ±0.14 


3.87±0.50 


2.08833 


1X30 




OSlRlS/10.4mGTC 


23.39±0.12 


1.60±0.18 


21.15017 


10X60 




OSIRIS/10.4mGTC 


24.21 ±0.14 


0.75 ±0.10 


28.49818 


5X180 




GMOS/8mGemini 


24.81 ±0.13 


0.43 ±0.05 


39.10159 


4X120 




OSIRIS/10.4mGTC 


24.77 ±0.13 


0.45 ±0.05 


39.47981 


5X180 




GMOS/SmGemini 


25 .24 ±0.15 


0.29 ±0.04 


44.08258 


4x180 




OSIRIS/10.4mGTC 


> 24.7 


< 0.48 


~180 


32x200 




OSIRIS/1 0.4mGTC 


26.90 ±0.14 


0.063 ±0.008 


1.17359 


17x180 




1.23mCAHA 


22.18±0.35 


4.88±1.57 


61.96267 


20X120 




SC0RP10/6mBTA 


25.17±0.35 


0.31±0.10 


0.29516 


3X300 




CQUEAN/2.1mMcD 


22.72±0.18 


2.96 ±0.49 


10.09449 


9X900 




RAT/2.0mLT 


24.01 ±0.13 


0.90±0.11 


39.12164 


5X60 




OSIRIS/1 0.4mGTC 


24.36±0.17 


0.65 ±0.11 


39.46336 


5X180 




GMOS/8mGemini 


24.61 ±0.09 


0.52±0.04 


0.30384 


3X300 




CQUEAN/2.1niMcD 


22.65 ±0.34 


3.16±1.00 


39.09432 


6x60 




OSIRIS/10.4mGTC 


24.73 ±0.42 


0.47 ±0.1 8 


39.44619 


7X180 


z' 


GMOS/SmGemini 


24.77 ±0.25 


0.45±0.10 


0.30745 


3X300 


Y 


CQUEAN/2.1mMcD 


> 22.5 


< 3.63 


37.45092 


32x60 


J 


NIRI/8mGemini 


> 23.4 


< 1.58 


28.46873 


44x60 




NIRI/8mGemini 


24.48 ±0.35 


0.59±0.19 



3.3 CAHA 1.23m 

The 1 .23m telescope is located at the German-Spanish observatory of Calar Alto (CAHA) in Almeria, 
is equipped with an optical imaging camera. The optical counterpart was detected in the VRI bands 1 
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days after the GRB trigger. The 1.23m was also used to calibrate the object field in BVRI bands by observing the 
Landolt fields RU149D and SA98 on Dec. 26 and 27, 2010, under photometric conditions. 

3.4 LT2.0m 

The Liverpool telescope is a 2.0m fully robotic telescope located at the observatory of Roque de los Muchachos 
on La Palma. Observations were carried out with the imaging camera RATCAM. The optical counterpart was 
detected in one epoch in i'-band at 10.09 d after the GRB. 

3.5 OSIRIS/GTC 

We acquired imaging data using OSIRIS at the Gran Telescopio de Canarias (GTC), a 10.4m telescope located at 
the observatory of Roque de los Muchachos on La Palma, Canary Islands, Spain. 

The observations started in r' band ~2 days after the burst, exposing for 30 s. A second r' observation was 
carried out ^21 days after the ganmia-ray event, where 5 exposures of 180 s were obtained. We furthermore 
obtained a late-time SED at 39 days in g' , r\ i' and z' bands and a last image at ~44 days in the r' band. 

The data of our two late epochs (at ~39 and ^44 days) were obtained at a considerable airmass (1.73 - 2.14) 
since the object was setting quickly after evening twilight. The data at ^21 and ^44 days were acquired with the 
Moon at ~54 and ~36 degrees, respectively, with an illumination of 83% and 21%, respectively. The SED at 39 
days was constructed in dark time. The observing conditions were good in all four GTC epochs. 

3.6 Gemini-North: NIRI and GMOS-N 

Late-time imaging of the optical counterpart of GRB 10 1225 A was conducted with the Gemini-North observatory 
on Mauna Kea/Big Island, Hawaii, on several occasions. On the night of Jan. 23, 201 1 we observed the field with 
the Near InfraRed Imager (NIRI) in the K' filter for 44 x 60 s exposures (2 x 30 s co-adds) before switching to 
the Gemini Multi- Object Spectrograph (GMOS-N) for 5 x 180 s exposures each in the r' filter. On the night of 
Feb. 1, 2011 we re-observed the field with NIRI in the J band for 32 x 60 s exposures (co-added), and finally on 
the night of Feb. 3, 201 1 we imaged the field in all four GMOS broad-band filters {g'r'i' z'). Since the source was 
setting, all exposures were taken at moderate to high airmass (1.5 — 2.5), although under relatively good seeing 
conditions. 
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3.7 BTA 6m 

A final late image was obtained using SCORPIO on the 6.0m BTA telescope, located at the Special Astrophysical 
Observatory, in Russia. The observation consisted of 20 x 120 s exposures using an / filter obtained on Feb. 25, 
2011 (60 days post burst) under good weather conditions and a seeing of 1.3 — 2.0 arcsec. 

3.8 Photometry of ground-based data 

The photometry of V, R and / data was done in a consistent way using a set of 15 comparison stars in the field 
of GRB 101225A calibrated with the Landolt fields taken on Dec. 26 and 27, 2011 by the 1.23m CAHA. For /, 
f and z' photometric calibrations was done with observations form the 2.1m Otto-Struve telescope at McDonald 
Observatory on Dec. 26, 2010, using the standard star Feige 34. Finally, photometry was derived from the 
other reference magnitudes using numerical transformations |[33ll . The magnitudes of the comparison stars in the 
different filters used for the optical observations are listed in Tab. [3] 

We performed aperture photometry using PHOT within IRAF taking an aperture radius equal to the Full Width 
at Half Maximum (FWHM) of the stellar point sources. In cases where the contamination by neighboring sources 
was not negligible, we did PSF photometry within IRAF. In Tab. [2] we list the final photometry for all UV, optical 
and IR data. 



Table 3: Magnitudes of calibration stars used for the optical photometry. A finding chart indicating the position of 
each reference star is given in Fig.|5] y,7?and/ are given in Vega system, the other bands are in AB. 



ID 


R.A. (J2000) 


Dec. (J2000) 


V 


R 


/ 


g' 


r 


•/ 

I 


z 


1 


00:00:48.13 


+44:38:18.9 


15.282 d= 0.009 


14.916 ± 0.049 


14.553 ± 0.052 










2 


00:00:52.02 


+44:37:46.6 


18.065 ± 0.097 


17.268 ± 0.055 


16.288 ±0.110 


18.62±0.23 


17.679±0.011 


17.121±0.010 


16.870±0.010 


3 


00:00:41.58 


+44:37:18.0 


16.794 ± 0.050 


16.421 ±0.051 


15.972 ± 0.062 


17.14±0.12 








4 


00:00:52.13 


+44:37:07.9 


18.789 ± 0.079 


18.335 ±0.133 


18.051 ±0.162 


19.26±0.19 


18.464±0.012 


18.325±0.011 


18.263±0.013 


5 


00:00:45.46 


+44:37:06.6 


14.995 ± 0.008 


14.633 ± 0.049 


14.235 ±0.051 










6 


00:00:46.18 


+44:36:52.4 


15.174 ±0.015 


14.838 ± 0.049 


14.456 ±0.051 










7 


00:00:55.37 


+44:36:36.8 


19.165 ±0.180 


18.268 ± 0.090 


17.356 ±0.188 




18.569±0.012 


17.905±0.011 


17.597±0.011 


8 


00:01:03.76 


+44:36:26.8 


16.594 ± 0.012 


16.161 ±0.049 


15.691 ± 0.090 










9 


00:00:48.48 


+44:36:19.3 


18.796 ± 0.059 


17.924 ± 0.069 


16.694 ± 0.092 


19.08±0.14 


18.602±0.012 


17.819±0.011 


17.516±0.011 


10 


00:00:47.98 


+44:35:57.8 


18.682 ±0.132 


17.901 ±0.107 


17.059 ±0.182 


19.41±0.31 


18.174±0.011 


17.645±0.010 


17.427±0.012 


11 


00:00:50.58 


+44:35:43.5 


18.505 ±0.218 


17.877 ± 0.060 


17.270 ±0.163 




18.088±0.011 


17.715±0.010 


17.537±0.012 


12 


00:00:51.59 


+44:35:19.1 


19.253 ±0.130 


18.893 ±0.131 


18.460 ± 0.273 


19.54±0.31 


19.059±0.014 


18.910±0.013 


18.852±0.020 


13 


00:00:43.29 


+44:35:13.1 


18.050 ± 0.066 


17.749 ± 0.085 


17.318 ±0.108 


18.34±0.16 


17.853±0.011 


17.729±0.010 


17.671±0.012 


14 


00:01:03.09 


+44:35:04.8 


16.828 ± 0.015 


16.321 ± 0.052 


15.826 ± 0.087 










15 


00:00:54.83 


+44:35:01.9 


17.516 ±0.032 


17.024 ± 0.066 


16.321 ± 0.082 


18.04±0.08 


17.148±0.010 


16.960±0.010 


16.853±0.011 
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Figure 5: Field with the secondary standards used for the photometric cahbration. The position of the optical 
counterpart of GRB 10 1225 A is indicated with an arrow. The photometric magnitudes of each of the reference 
stars are given in Table |3] The field of view is 6 arcmin x 4 arcmin. 

3.9 Preimaging 

Preimaging of the field was available from the archive of the 3.5m Canada-France-Hawaii Telescope (CFHT), 
obtained with the MegaPrime/MegaCam for the Pan- Andromeda Archaeological Survey (PAndAS, 1341 ). We 
combined 3 x 500 s exposures obtained under very good conditions in g' and i' bands and derive 3cr limiting 
magnitudes for these exposures of i' > 25.5 and g' > 26.9. At the position of GRB 10 1225 A, we tentatively 
detect an object with g' = 27.2 ± 0.5 (2cr, see Fig. [6]) which is consistent with our late host detection (see Sect. 

3.10 Host detection 

At three epochs, June 9, June 28 and 29, we obtained deep imaging of the field of GRB 101225A in g' and r' 
with OSIRIS at the 10.4m GTC under good conditions and seeing around 1.0 arcsec. The total exposure time is 
8400 s in g' and 6400 s in r' band. We detect an unextended object at the position of the GRB in both bands at 
g' = 27.21 ± 0.27 (-3cr) and r' = 26.90 ± 0.14 (-7a). These values are clearly above the extrapolation of the SN 
light curve to 180 days and we therefore propose this object to be the host galaxy of GRB 101 225 A. At a redshift 



20 




Figure 6: Pre-imaging exposure in g' band obtained with the 3.5m CFHT. The field of view is 60 arcsec x 
40 arcsec, North is to the top and East to the left of the image. The circle indicates the position of the optical 
counterpart of GRB 101 225 A, where we see a 2cr detection of what could be the host galaxy. 

of z = 0.33, the absolute luminosity is Mabs,g = —13.7 mag, ~ 2mag fainter than the faintest GRB host detected 
so far (XRF 060218 [12]) and one of the faintest galaxies ever detected at that redshift. The g' — r' color of the 
host is consistent with what we would expect for a late-type star-forming galaxy, in line with what is seen for other 
long GRBs, although we note that the large errors allow other galaxy types as well. The blue color would also 
argue against the source being a late detection of a cool component from the event itself (see [3]). 




Figure 7: Imaging of the host using OSIRIS/GTC 6 months after the GRB. Panel a shows the stacked image 
in r' , panel b the same in the g' band. The field of view is 50 arcsec x 65 arcsec. The circle indicates the position 
of the host galaxy. 
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4 Optical spectroscopy 



We obtained a spectrum of the optical counterpart 51 h after the event using OSIRIS on the 10.4m GTC telescope 
on La Palma (Spain). Two spectra of 1800 s exposure time each were taken with grism 300B (R=325, wavelength 
range: 3500 — 7000 A) under moderate to high airmass (1.26 and 2.05, respectively). The spectra were reduced 
and combined with standard tasks in IRAF and flux-calibrated with the spectrophotometric standard G191-2B2 
observed in the same night. The continuum is clearly detected, but the spectrum shows no obvious absorption or 
emission lines. 

We performed a search for any possible redshift solution in the range z=0.1 - 0.6. For this, we stacked the 
normalized spectrum at the position of the strongest emission lines [Oil] AA 3727, 3729 A, [OIII] AA 4959, 5007 
A, H/3 and Ha at redshifts between 0.1 and 0.6 in steps of 0.005. We find no indication of any emission feature in 
the stacked spectrum at any redshift in this interval. 

The limits on the detection of Ha [OIII] and [Oil] emission from the host galaxy are < 5 x 10~^^ erg cm~^ s~^, 
< 2.3x 10~^^ ergcm"^ s~^ and < 3x 10~^^ ergcm"^ s~^ (3a) respectively. We can also put a limit on the de- 
tection of Ha at z = of < 2x 10~^^ ergcm"^ s~^. The flux-calibrated spectrum with the position of typical 
emission lines from the host shifted to a redshift ofz = 0.33 is shown in Fig. [5] 




4000 5000 6000 7000 8000 9000 10000 

observed wavelength in A 

Figure 8: Flux-calibrated spectrum of the optical transient from GTC/OSIRIS 2.1 days after the GRB. The 

error spectrum is plotted in blue. The red lines indicate the position of normally strong emission lines from the 
interstellar medium at a redshift of z = 0.33, none of the lines is detected in our spectrum. 

On the night of 201 1 Feb. 04 we observed the optical counterpart with the Low-Resolution Imaging Spectrom- 
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eter 1351 on the Keck I telescope during local twilight. Two undithered observations of 600 seconds each were 
acquired using the 1.0 arcsec slit and the D500 dichroic at a position angle of 86.5 degrees. On the red arm we 
used the 600/7500 grating and binned the CCD along the spatial direction (2x1 binning); on the blue arm we used 
the 600/4000 grism and binned the data along both spatial and spectral axes (2x2 binning). Due to twilight there is 
no evidence of a trace in the blue spectrum (and no source is detected in GMOS ^'-band imaging from the previous 
night). A faint continuum trace at the expected position of the transient is identified on the red side in the second 
(less twilight-affected) exposure from 7160 to 8000 A with no absorption or emission features visible. 

5 Modeling the UV to NIR spectral energy distribution 

The early evolution of the UV/optical/IR (UVOIR) counterpart is very unusual for a GRB afterglow. Instead of a 
power-law spectrum with a negative spectral slope, it had a very blue counterpart, following what seemed to be a 
power-law with a positive spectral slope |36|. Furthermore, the counterpart stayed bright during the first days and 
then decayed (31] with a strong color change, transforming into a very red counterpart two weeks after the trigger 
1381 . We interpret this early evolution as being produced by the expansion and cooling of a blackbody (BB), as 
shown in Section |5J1 

The simple BB evolution is not valid any more for the emission beyond ~ 20 days after the trigger. At that 
time we observe a flattening of the light curve, while the very red color is preserved. This late evolution can be 



well-described with the presence of a supernova component, as described in Section 5.2 

5.1 Early time evolution 

For the modeling of the UVOIR spectral energy distribution (SED), we use the photometry presented in Tab. [2] 
together with some of the data points extracted from the literature |[33[37l[38l|39l|40l|4Tl|42|. All magnitudes are 
corrected for a Galactic extinction of Ay = 0.33 mag and transformed from magnitudes to flux densities. With 
this data set we are able to derive a set of 12 SEDs ranging from 0.07 to 40 days after the trigger. 

The early optical SEDs are well fitted by using an expanding and cooling blackbody of the following form (in 
frequency space): 

F.{Jy)-10 ^-j (1) 
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Here, the factor 10^^ is used to convert W/m? /Hz to Jy. R is the radius of the emitting black body (which we 
assume to be spherical), D is the luminosity distance to the object, z the redshift and Tohs is the observed blackbody 
temperature (the rest-frame temperature would be Trest = ^o6s(l + ^))- The other physical constants are: c is the 
speed of light, h Planck's constant and kB Boltzmann's constant. For simplicity we assume a blackbody with an 
emissivity of 1. The blackbody succeeds in reproducing the data up to 10 days, without any intrinsic extinction or 
additional emission component, after which another component becomes dominant (see Section [5^ . 

Table 4: Measured values for the blackbody evolution. Values in brackets are extrapolated from the BB evolution 
due to a limited amount of data points in those SEDs. 



Epoch 


Observed temperature 


Normalisation constant 


(days) 


(K) 


(l + z)4^1026(g)' 


0.07 


(43 000±8 000 


1.7±1.5) 


0.17 


40 000±6 000 


1.8±1.2 


0.3 


35 000±3 950 


2.2±0.8 


0.6 


25 340±5 440 


4.6=b2.5 


1.1 


20 900±1 770 


5.2±1.2 


2.0 


15 000±1 090 


8.9±1.8 


3.0 


14 260±1 760 


8.00±3.4 


5.0 


(11 300±2 000 


10.0±4.0) 


10.0 


(6 000±2 000 


14.0±4.0) 


18.0 


(5 000±1 000 


18.0±10.0) 



From the fits to the SED evolution and allowing for a second-order fit, we get the following evolution of the 
normalisation constant: 

log ^10^^7r(l + zf j = (0.70 ± 0.04) + (0.46 ± Om)log{t) - (0.01 ± Om)loglO{tf (2) 

where t is the time in days. The temperature evolution (in K) can be described by: 

log {Tabs) = (4.342 ± 0.017) - (0.395 ± 0.016)%(t) - (0.11 ± 0.02)log{tf (3) 

Figure [9] shows the temporal evolution of the normalisation constant and the temperature. The normalisation 
can be sufficiently described by a linear evolution in log-log space and therefore the second order term in eq. (2) 
can be neglected. For the temperature, we need an additional second order term to obtain a reasonable fit to the 
data. The temperature and normalization evolution fits very well to our theoretical model of the afterglow (see Sect. 
[8] and Fig. [19]). Once the redshift is known the normalisation constant can be transformed into physical values in 
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the rest- frame of the object which is shown in Fig. 
Mpc using a ACDM cosmology with Hq = 71, 



10 For this we assume z 
= 0.27 and = 0.73. 



= 0.33 (see Section [5^, or 1661.1 
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Figure 9: Temporal evolution of the normalization constant and temperature during the UVOIR BB phase. 

The evolution of the normalization constant can be described by a linear evolution in log-log space while the 
temperature evolution requires a second order term. 



5.2 Late evolution and SN template fitting 

The late evolution of the light curve requires a component in addition to the evolving BB described in Section 



5.1 A pure BB fit to the late data gives a bad fit with x /d.o.f. = 10.45/5 and is therefore ruled out with 
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Figure 10: Temporal evolution of the physical radius and temperature of X-ray and UVOIR BB. This figure 
is identical to[9]but showing the evolution of both BB components together. 

a 94% probability. We assume the late SED to have an additional component from a SN and fit the SED with 
several SN templates. These fits are also used to estimate the redshift of GRB 10 1225 A, which we could not obtain 
spectroscopically (Section]?]). 

To determine the redshift of GRB 101 225 A we use the SED at 40 days after the burst where the contamination 
from the BB is negligible and where we have detections in 7 different bands. Given the steep slope in the blue part 
of the spectrum, we convolve the response of each filter with the spectral shape of the templates. This is particularly 
important for the r'-band observations performed from GTC and Gemini at a very similar epoch, which show a 
significant difference in flux density. The filter of GTC reaches slightly redder wavelengths, and the difference in 



flux densities can be well explained by a very steep slope due to a SN feature as shown in Fig. 12 



We obtain templates for different core-collapse supernovae from the literature We exclude SN la from the 



http://supernova.lbl.gov/^nugent/nugent-templates.html^ 
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0.1 1.0 10.0 

observed time (days) 

Figure 1 1 : Evolution of the velocity of the UVOIR black body according to the result of our modeling. The 

plot shows the evolution at the best-fit redshift ofz = 0.331q q4 (red dashed line and grey dotted lines) as well as 
three other redshifts. At z = 0.33, the velocities derived during the BB phase range between 0.25 — 0.005c, typical 
of accelerated material during a supernova explosion. 

analysis, as we do not expect high-energy emission such as detected for GRB 10 1225 A for those events. The 
template for each SN prototype was interpolated to the time of the SED for a range of redshifts (see references 
in Table [5]). In the particular case of SN 1998S the templates were created by combining ground-based and HST 
spectra. 

In order to obtain the best-fit redshift, we also need to consider the time evolution and maximum of the SN 
light curve which is expressed as the stretching factor (a stretching factor of s = 1 corresponds to a time evolution 
identical to the corresponding template SN). To this end, we iteratively fit the SED and the light curve starting with 



5 = 1 and a first fit to the SED in 7 bands (see Fig. [12]). With the best-fit z, a SN light curve is derived and fitted 
to the real light curve which gives a new stretching factor. This process continues until the solutions for s and z 
converge. We also allow for a scaling of the flux of the SN during the fit. The best fits for s and z for each SN 
template are displayed in Table [5] 

The absolute best fit is obtained with a SN 1998bw template, a broad-lined Type Ic that is the canonical template 
for GRB -related supernovae. For this case we obtain a redshift of ^ = 0.33to;o4 and 5 = 1.25 ± 0.15. The other 
core-collapse SNe we tested (excluding the Type II SN 1998S, which clearly does not fit our SED) give redshifts 
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Table 5: Fits of the SED with SN templates. A pure BB fit (not listed in this table) is rejected with a probability of 
94%. 



SN Type 


SN template 


best fit z 


stretching factor 




rejection prob. 


reference 


Ib/c 


1999ex 


'-'•■^^-0.10 


1.1±0.15 


5.80/5 


67% 


na 


Ic broad-lined 


1998bw 


33+^-°'^ 

V.JJ_0 04 


1.25±0.15 


3.60/5 


39% 




II 


1998S 


'-'•-''-'-0.08 


1.0 ±0.15 


44.5/5 


100% 


isiiiTiiiii 


UN 


200 laj 






> 50.0/5 


100% 


im 


IIL 


1985P 


Q .0+0.04 


1.3 ±0.15 


6.1/5 


70% 


eqi 


IIP 


1999em 


41+0-05 

'-'•^^-0.03 


1.0 ±0.15 


6.3/5 


72% 


f5T1 



between z = 0.31 and z = 0.50. We therefore use z = 0.33 as reference in this work. Figure 12 shows the fitted 
SED with the template at z = 0.33to.o4 Figure 13 the best fits for all other SN templates as well as a simple 
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Figure 12: Fit of the day 40 SED to a SN 1998b w template. The observations are in black and the best fit, with 
a redshift of z = 0.33, in red. The gray lines represent the template at the different redshifts (steps of 0.01) within 
errors. 



5.3 Luminosity and stretching factor of the SN associated with GRB 101225A in context 
of other GRB-SNe 



We undertake a more general comparison to SNe associated with GRBs by following the formalism of f52l. 
These authors used a SN 1998bw template light curve to fit late bumps in GRB afterglow light curves, modifying 
the template by increasing or decreasing the luminosity at peak (the parameter /c, with /c = 1 implying a peak 
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Figure 13: Fit of the SED at 40 days to SN templates other than a SN 1998b w-hke (as listed in Tab. [5| The fit 
shown is the one with the best-fit stretching factor for each template. We also plot a simple BB with 3650 K, the 
temperature at which we obtain the best fit. 

luminosity identical to that of SN 1998bw), and stretching or compressing the light curve in time while retaining 
the overall shape (the parameter s, again, 5 = 1 implies the temporal evolution is identical to that of SN 1998bw 
in the same band). This procedure also included the creation of synthetic templates by interpolating between the 
SN 1998bw light curve in different filters, and taking into account the cosmological K-correction. Nearly all 
GRB-SNe were well-fit by the SN 1998bw light curve template. For GRB 101225A we have s = 1.25 ± 0.15 and 
k = 0.08 ± 0.03 according to the designation of 15211 using the SN 1998bw light curve at z = 0.33. 

Ferrero et al. | 53 1 analyzed SN 2006aj associated with XRF 060218, and placed it into the /c — s context. They 
employed the line-of-sight extinction values derived by 1541 to derive intrinsic k values. To place the SN associated 
with GRB 101 225 A into the /c — 5 context, we fit the light curve analogous to 1521 , and use the sample of l53l as 
well as additional events as a comparison. The complete data are presented in Table [6] 

GRB 990712 has been analysed again with additional data. We find no evidence for host extinction. For GRB 
021211, a re-analysis of the afterglow SED finds no evidence for host extinction, the value from l53l thus remains 
unchanged but now counts as extinction-corrected. For GRB 040924, we use the k and s values from [9] and 
correct k with the extinction found by |54|. For GRB 050525 A, we use the uncorrected k value from |53|, and 
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correct it with the extinction found by [11]. For XRF 050824, we use the uncorrected k value from |55|, and 
correct it with the extinction found by [11]. GRB 060729 is analysed in (561. GRB 080319B is analysed in |57 |. 
GRB 090618 has been analysed for this work, using the data set of [18]. We were not able to derive a good SED 
for this afterglow, therefore the k value has not been corrected. 



As can be seen in Fig. 14 the SN associated with GRB 101225A is significantly fainter than any other known 
GRB SN (with the SN associated with GRB 040924 being the most similar, but this event is only marginally 
detected). At the same time its temporal evolution is similar but slightly slower than most known GRB-SNe, 
though not by a large amount. It is also fainter than two well-studied Type Ic SNe, SN 19941 and SN 2002ap (see 
1531 for discussion), the latter being broad-lined, but not associated with a GRB. 
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Figure 14: Luminosity factor k and stretching factor s of SNe associated with GRBs. Filled symbols have 
been corrected for host-galaxy line-of-sight extinction, non-filled symbols have not. We label several well-studied 
nearby GRB-SNe, as well as two "canonical" Type Ic SNe, SN 19941 (Ic) and SN 2002ap (broad-lined Ic unasso- 
ciated with a GRB). The GRB 101225A SN is fainter than all these events. 



6 Discussion on the redshift of GRB 101225A 

Determining the distance scale at which GRB 10 1225 A occurred is crucial to understand the energetics and get a 
clear picture of the physics involved in this event. In this section, we present several independent arguments to 
strengthen our redshift estimation. 

The first strong limit on the redshift comes from the UVOT detection in uvw2 which implies a redshift lower 
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Table 6: Luminosity Factor k and Stretching Factor s for GRB SNe 



LrKr> 


K 


s 




Comment 


Reference 


y /{jZZo 


n /i 1 j_ n on 
U.4i ± U.zy 


1 /I IT in 
i.4o ± U. 


nQ 
.yo 


uncorrected 


1 


QQn'71 

yyu / iZ 


U.oU ± U.iU 


U. / ± U. 


.U ( 




This Work 




r, or+0.44 


1 /in -k n 






ID J II 


ni 0091 


1 oi::+2-82 
i.OO_o.79 


n f\Q -u n 
u.oy =t u. 


91^ 
.ZD 




iR^ 

ID J II 


Ui i izi 




U.oU ± U. 


no 
.Uz 




ID 3 II 


020405 


n nn+O.lS 

0.90^0.11 


0.97 ± 0. 


.07 




1531 


\JZ\Jy\Jj 


u.oz ± u.uy 


n no -U n 
U.yz ± U. 


.Uo 


uncorrected 




Uzizi i 


U.4o ± U.zi 


n nn j_ n 

u.yy ± U. 


1^ 
.zo 




tjjj , Inis WorK 


030329 




0.85 ± 0, 


.10 




ED 


031203 


-1 90+0. 18 


1.09 ±0, 


.07 




EH 


040924 




1.37 ±0, 


.97 




[10], EH 


041006 




1.38 ±0, 


.06 




ED 


050525A 




0.77 ±0, 


.04 




f53|,[ll] 


050824 


-'-•^'-'-0.56 


0.52 ±0, 


.14 




155 1, [11] 


060218 


0.74 ±0.01 


0.69 ±0, 


.01 




153J 


060729 


1.04 ±0.03 


0.86 ±0, 


.02 




m 


0803 19B 


^•'^'-'-0.90 


0.89 ±0, 


.10 




m 


090618 


0.58 ±0.05 


0.99 ±0, 


.05 


uncorrected 


This Work, [19] 


101225A 


0.08 ±0.03 


1.25 ±0, 


.15 




This Work 



than z = lA 1581 . As an independent test, the analysis of the absorption in the X-ray spectra (see Sect. 2 and 
in particular Fig. |4]) imposes an upper limit on the redshift of 0.5 within 99% confidence and 0.35 within 90% 
confidence. 

The most restrictive redshift estimate comes from the SN fitting at 40 days after the burst (see Sect. 5.2), 
where we find a best fit for a broad-lined SN Ic, such as SN 1998bw, the prototype of GRB-related SNe, at 
z = 0.33^0 04. For the rest of the SN types that give a reasonable fit, we obtain similar values, always resulting in 
a redshift between 0.21 and 0.50. Independently of the SN fit, we can compare the brightness of the bump in the 
fight curve with the dimmest known SN (SN2008ha, which peaked at Mr = -14.5 ± 0.3, |59 |) and take that as a 
lower limit on the redshift. In this case, we would expect GRB 101 225 A to be located at a redshift larger than 0.1. 

From the SED fit of the first days, we know that the evolution is well-described by a simple BB. Depending on 
the distance at which the object is found, we can derive different radii and expansion velocities. For an explosion 
of this type, we expect expansion velocities larger than ~ 10^ km s~^, which would be at the edge of a stellar 
wind regime and, if similar to a SN explosion, of the order of 10^ km s~^. We cannot, in principle, rule out 
higher velocities of the ejecta. Under the assumption that the ejecta should not be traveling at velocities larger 
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than 100,000 km (0.3c) and lower than 1,000 km (0.03c), we can estimate a range of reasonable velocities 
between z ~ 0.20 and z ~ 0.60 (see Fig.[TT]). We note that there might be additional effects such as variations in 
the transparency of the BB that could introduce variations in this simplified analysis. In any case, the evolution of 
the BB would be hard to get by any source within the Local Group. At a redshift of 2; = 0.33, the velocity of the 
blackbody would have evolved from ~ 70, 000 km s~^ at the time of our first SED to ~ 2, 000 km s~^ at the time 
of our last blackbody-dominated epoch, nicely matching the requirements. The fact that the blackbody detected for 
XRF 060218/SN 2006aj was very similar to the one found for GRB 101225A when placed at a redshift of z = 0.33 
(see Section [7]) adds additional evidence for the validity of the redshift estimate. 

Concluding, the proposed redshift for GRB 101225A of z = 0.33to.o4 supported by several independent 
arguments, and can be considered as a firm reference when studying the physical processes involved in the event. 
We do not find any evidence that would indicate a redshift smaller than z = 0.2 or higher than z = 0.5. 

7 Comparison between GRB 101225A and other GRBs with SNe and BB 
components 

In Table |7j we compare several properties of GRBs associated with SNe and without a "classical" afterglow com- 
ponent. All of them are subluminous compared to the average long-duration GRB with Eigo around 10^^ — 10^^ 
erg. GRB 101225A lies on the lower end of the energy output from normal long-duration GRBs. Among those 
nearby GRB -SNe without a classical afterglow, there is a class of very long duration GRBs with very low Epeafe 
values, all of them showing a thermal component in X-rays. XRF 060218 and XRO 080109 also had a thermal 
component at optical wavelengths during the first few days [13, 22]. For XRF 1003 16D, no optical counterpart 
was detected before the onset of the actual SN due to high intrinsic extinction in the host galaxy. 

XRF 060218 [13] shows a similar early behavior to GRB 10 1225 A. We compare the early UVOT light curve of 
XRF 060218, which we obtained from the UVOT catalogue (6OI, to the fight curve from GRB 101225A by shifting 
that of XRF 060218 to z = 0.3, including a i^-correction. To obtain the i^-correction, we use XSPEC assuming 
a blackbody spectrum with kT ~ 3.7 eV. The temperature was determined from the best-fit model of a BB to a 
SED of XRF 060218 taken at 120ks after the trigger [13]. Using this BB spectrum, we determined the expected 
flux density in the observed frame for each filter and at z = 0.33. The ratio of these two flux densities was taken to 
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Table 7: GRBs with SNe but without afterglows 



GRB 


z 


TC? 


Tgo 




F 


HR(50-100) 


Radio? 


SN My 


Host Mb 








(s) 


(keV) 


(erg) 




(mag) 


(mag) 




980425 


0.0085 


No 


23.3 


55±21 


8.1 xlO"^^ 




Yes 


-19.42 


-17.6 


031203 


0.105 


No 


30 


158±51 


3x10^^9 




Yes 


-20.39 


-21.0 


060218 


0.0331 


Yes 


-2100 


4.9 


6.2x10"^^ 


0.835 


Yes 


-18.76 


-15.9 


080109* 


0.0065 


Yes 


-400 


low 


2x10^^6 




Yes 


-16.7 


-20.7 


1003 16D 


0.059 


Yes 


>1300 




3.1 xlO"^^ 


0.891 


No 


-18.62 


-18.8 


101225A 


0.33 


Yes 


>2000 


38±20 


> 1.4 X 10^^ 


1.06 


No 


-16.9 


-13.7 



*No 7-rays observed, numbers derived from X-rays. 

TC? refers to the early thermal component, mostly attributed to a supernova breakout. 
HR is the hardness ratio, defined as the ratio of channels (50 — 100 kev)/(25 — 50 keV) 
SN My is the SN peak absolute magnitude in V 
Host M^ is the host absolute magnitude in B 



be the i^T-correction for the specific filter, which is — 2.20 for all filters. Figure \T5\ compares the flux-density light 
curves of the 3 UV filters of both GRBs. For both GRBs the light curves were corrected for Galactic extinction. 




Figure 15: UV light curves in flux density for GRB 101225 A and XRF 060218. The fight curves for XRF 
060218 have been shifted to z = 0.33 for direct comparison. The colored solid shapes connected by solid lines are 
for GRB 10 1225 A, open symbols and dotted lines for XRF 060218. Circles are uvwl, squares are uvm2, triangles 
are uvw2. 



We furthermore performed a fit of the early UVOIR SEDs of XRF 060218 and SN 2008D, in a similar way to 
what we did for GRB 101225 A. XRF 0602 18 also seems to follow a blackbody evolution at early times. However, 
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the SN starts to dominate already around 3 days after the burst (see Fig. [16]), Hmiting how long the study of the 



evolution is possible. As can be seen from Fig. 17 the evolution of the blackbody is not very different from what 
we see in GRB 10 1225 A, although the radius expansion is slightly steeper. SN 2008D does not have a measurable 
thermal component in the early X-ray data. In the optical, the cooling BB dominates until about 4 days after the 



event when the onset of the SN was observed (see Fig. 16). The temperature of the UVOIR BB of SN 2008D is 



lower (- 30 000 K) than those of XRF 060218 and GRB 101225A. The radius evolution is considerably steeper 



than for those two events (Fig. 17 ). The UVOIR BB emission of SN 2008D is therefore likely due to the cooling 
of the initial shock breakout. 

8 Temperature evolution and geometry of the UVOIR BB 

The interaction of an ultra-relativistic, well-collimated jet (having a half-opening angle Oj^o ~ 10°) with the 
previously ejected broad torus can explain the power-law and thermal component in the X-rays of GRB 10 1225 A. 



In Fig. 18 we show four different stages in the evolution of the jet and its interaction with the CE shell. 

Taking the results from CE- simulations, we assume that the density of the common envelope ejecta in a narrow 
funnel around the symmetry axis of the system is much lower than elsewhere in the CE-ejecta. This funnel has 
an opening half-angle 6f ~ 2°, which permits the passage of the 7-radiation generated by internal shocks in the 
ultra-relativistic jet. Since Of < Oj^o, most of the jet beam hits the inner boundary of the CE-ejecta at a distance of 



RcE,in — 2.5 X 10^^ cm (panel a of Fig. 18 ), while only a small fraction of the beam propagates through the funnel 
until it reaches the outer radial boundary of the CE-ejecta (RcE,out — 2.1 x 10^^ cm, panel b). As the central 
jet spine progresses through the ejecta funnel, it interacts with the lateral walls, giving rise to mass entrainment in 
the jet. The additional baryon load and the shear with the funnel walls decelerates the ultra-relativistic jet spine 
extremely quickly and, therefore, no regular afterglow signature is produced. Furthermore, most of the jet beam 
hits regions of the CE-ejecta which are much denser than the central funnel and, thus, its propagation across the 
ejecta is much slower than that of the central spine. As the jet decelerates to moderately relativistic speeds it 
expands laterally. 

The X-ray emission is attributed to the shocks (forward and reverse) produced as the ultra-relativistic jet im- 
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Figure 16: SED fit of the early UVOIR data of XRF 060218/SN 2006aj (top) and SN 2008D (bottom). The last 
epoch in both plots already begin to show a strong contribution of the SN. These figures can be directly compared 
to Fig. 1 of the main paper. 

pinges against the inner radial boundary of the CE-ejecta. In our model the cross-sectional radius of the CE-ejecta 
funnel (Rf = RcE,in sin Of ~ 10^^ cm) sets the almost-constant size of the observed X-ray source (blue dots in 



Fig. [TO]). 

An approximate evolution of the jet spine as it crosses the ejecta funnel can be obtained by applying the model 
of 1631 , which is based on solving a system of four ordinary differential equations (as a function of the observer's 
time t) for the Lorentz factor, the rest-mass, the radius and the jet half-opening angle, respectively: 
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Figure 17: Temporal evolution of the radius and temperature for the pre-SN epochs of XRF 060218/SN 
2006aj (top) and SN 2008D (bottom). The dotted line shows a polynomial fit to the temperature and radius 
evolution. This figure can be directly compared to Fig. [TO 
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dT dm 

— = (4) 

dt dt Mjet + em + 2(1 - e)Vm 

dm dR , 

-— = -—27ri?^(l -cos (9)n(i?) (5) 

LLL LLL 

dR 

~dt 



PcT{r + - 1) (6) 



d^ _ c,(r + vT^^) 

dt R ' 

where n(i?) is the particle number density, j3 = Vl — F^^ is the velocity, e is the shock-generated thermal 



energy (in the comoving frame) that is radiated, and Cg = \J (4r + l)(r^ — l)(3r(4r^ — 1))"-^ is the sound 
speed. The temperature of the jet spine is derived by assuming that the pressure is dominated by the photon 
radiation, so that P(t) = aT^(t)/3, with the radiation constant a = 7.56 x 10~^^ ergcm~^K~^, and noting that 
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P = (r^ — 1) / {3T)nmpC^ , (rrip being the proton mass) according to 1631 . 

For the jet spine propagating through the CE-ejecta funnel we may take the following initial conditions, corre- 
sponding to an ultra-relativistic and cold jet: Tin ^ 100, M^^ 8 x 10~^Mq, Rin = RcE,in = 2.5 x 10^^ cm, 
and 6in = f = 2° . A moderate radiative efficiency e ^ 0.4 is assumed. The initial jet mass is estimated from 
the observed lower limit of E^^iso and the initial jet half-opening angle as Mjet = Ejet/{{iX ~ 1)(1 ~ ^) + 
where Ejet = Ej^iso{l — cos 6f) 8.5 x 10^ erg. We assume that the CE-ejecta have a declining particle number 
density of the form n{R) = nin{RcE,in/ ^ with riin = 4.9 x lO^^cm"^. After t 0.05 days the jet spine 
reaches RcE.out at a velocity of v jet, out ^ 0.25c, an opening half-angle of jet, out — 70°, and having plowed 
1.5 X lO-^M© of the CE-ejecta. 



The temperature and radius evolution of the UVOIR BB (panel c in Fig. 18 ) can be modeled as the result of 
the further deceleration of the sub-relativistic jet emerging trough RcE^out- Again, we use equations [4j|7] but now 
considering a uniform medium with a particle density of riext — 1.6 x 10^ cm~^ and a higher radiative efficiency 
of e = 1. As initial conditions we take the terminal values of the velocity, mass, radius and opening half-angle of 



the previous evolution through the CE-ejecta funnel. The result can be seen in Fig. [19] (after t ~ 0.05 days). Small 
variations of the initial variables can also roughly fit both the radius and the temperature evolution. We note that, 
even though the theoretical model does not predict an exact power-law for the radius evolution, the deviations from 
the observed R{t) ex t^-^^ are rather small. We also point out that the temperature evolution is not compatible with 
a single power-law as the observational data suggests. 

Finally, we speculate on the fate of the jet wings, i.e., the fraction of the jet that does not cross the CE-ejecta 
funnel. The energy contained in these wings, ~ 10^^ erg, is transferred to the CE-ejecta, and can break out 
through RcE,out on scales of ~few days, almost isotropically, but at relatively low temperatures compared with 
the emerging jet spine. Hence, its observational signature will be probably hidden by the emerging SN light curve 



(panel d) in Fig. 18 
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Figure 19: Radius (panel a) and temperature (panel b) evolution of a sub-relativistic jet. The initial evolution 
up to 0.05 days (grey shaded area), is computed with a modified version of the 1 63 1 method to incorporate a shell 
profile in which the number density decays with radius as n = nin{RcE,in/R)'^, with riin = 4.9 x 10^^ cm~^. For 
this part of the evolution, the initial data are: Rin = RcE,in = 2.5 x 10^^ cm, Ejet — 8.5 x 10^^ erg, F^^ ^ 100, 
Oin = 6>/ = 2°, and e:^ 0.4. 
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